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Abstract. We apply a heuristic method based on counting points over finite fields to the 
Poincare center problem. We show that this method gives the correct results for homoge- 
neous non linearities of degree 2 and 3. Also we obtain new evidence for Zoladek's conjecture 
about general degree 3 non linearities. 



Introduction 

In 1885 Poincare asked when the differential equation 

/ = x+p(x,y) = _ P(x,y) 

y + q(x,y) ' Q{x,y 

with convergent power series p(x, y) and q(x, y) starting with quadratic terms, has stable 
solutions in the neighborhood of the equilibrium solution (x, y) = (0,0). This means that in 
such a neighborhood the solutions of the equivalent plane autonomous system 

x = y + q(x,y) = Q(x,y) 
y = -x- p(x, y) = -P(x, y) 

are closed curves around (0,0). 

Poincare showed that one can iteratively find a formal power series F = x 2 + y 2 + f2(x, y) + 
H( x i V) + • • • such that 



( \ CO 



3=1 

with Sj polynomials in the coefficients of P and Q. If all sj vanish, and F is convergent then 
F is a constant of motion, i.e. its gradient field satisfies Pdx + Qdy = 0. Since F starts with 
x 2 + y 2 this shows that close to the origin all integral curves are closed and the system is 
stable. Therefore the s,-'s are called the focal values of Pdx + Qdy. 

Poincare also showed, that if an analytic constant of motion exists, the focal values must 
vanish. Later Frommer |Fro_34j proved that the systems above are stable if and only if all 
focal values vanish even without the assumption of convergence of F. (Frommer's proof 
contains a gap which can be closed |vW05j ) 

Unfortunately it is in general impossible to check this condition for a given differential equa- 
tion because there are infinitely many focal values. In the case where P and Q are polynomials 
of degree at most d, the Sj are polynomials in finitely many unknowns. Hilbert's Basis The- 
orem then implies that the ideal Iqq = (si,S2, •••) is finitely generated, i.e there exists an 
integer m := m{d) such that 

si = s 2 = ■ ■ ■ = s m ( rf ) = =J> Sj = Vj. 



Date: February 1, 2008. 



1 



2 



V. BOTHMER 



This shows that a finite criterion for stableness exists, but due to the indirect proof of Hilbert's 
Basis Theorem no value for m(d) is obtained. In fact even today only m(2) = 3 is known. 



Zoladek |ZoI95| and Christopher |Chr05bj showed that m(3) > 11. Our experimental data at 



this moment recovers m(3) > 10. 

The proof for m(2) = 3 is conceptually simple: Compute the first 3 focal values as polynomials 
in the coefficients of P and Q under the assumption deg(-P) = deg(Q) = 2. The 3 polynomials 
cut out an algebraic variety in the space of all differential equations of degree 2. Then 
decompose, by hand or by computer, this variety into its irreducible components. For each 
component prove that all its differential equations have a constant of motion. By works of 
Dulac Dul08 and Schlomiuk |Sch93| we know that in this case integrating factors of Darboux 
type, i.e of the form ]X =1 F" 1 with F{ polynomial always suffice to find constants of motion. 

For d = 3 this approach is not feasible because the polynomials Sj are very large. They 
involve 14 variables and are of weighted degree 2j. For example S5 has already 5348 terms 
and takes about 1.5 hours on a Powerbook G4 to calculate. The polynomials Sj, j > 6 can 
at the moment not be determined by computer algebra systems. Even if we would somehow 
obtain these polynomials, it is extremely difficult to decompose the resulting variety into 
irreducible components. Even I5 = (s\, . . . , 35) can not be decomposed by current systems. 
So for d = 3 only partial results are known. 

One new feature of the center problem in degree d = 3 is that a new type of differential 
equations with a center appear, namely the rationally reversible centers. Z olade k has cl assified 
all rational reversible cubic systems which are not of Darboux type in |Zol94j and [Zol96 . 
He also conjectured that all cubic systems with stable solutions near the origin are either 
of Darboux type or rationally reversible. This conjecture has been ver ified on several linear 
subspaces of the space of all differential forms, for example in |CRZ97| and |Chr05a| . In this 
paper we provide additional statistical evidence for Zoladek's conjecture in the whole space 
of Poincare differential equations up to codimension 7. 

Our main tool is a statistical method of Schreyer lvBS04| to estimate the number of compo- 
nents of the locus X{ where the first i focal values vanish. The basic idea is to reduce the 
equations Sk modulo a prime number p and count the number of F p -rational points of Xi. 
By the Weil Conjectures }Wei49] . which were proved by Delinge |Del74j . we know that the 
fraction of points 

number of F„ rational points on Xi 
l P {Xi) ■- 

is equal to 

T 



r y~j + higher order terms 

where c is the maximal codimension of a component of X and r is the number of components 
with this codimension. 

To evaluate the st at all possible points is not feasible, but by using a large number of random 
points one obtains an approximate value of J P (X) that can be used to estimate r and c. 

In an appendix to this paper which is joint work with Martin Cremer we show how one can use 
Frommer's algorithm to evaluate the s^'s in a given point over a finite field ¥ p without knowing 
the explicit polynomial equations for s&. We describe a C++ implementation following 



|Kro34j . |Mor00| and |H0l| that can evaluate 3.000.000 points per hour on a 450 MHz machine. 



As a slight improvement to Schreyer's method we also look at the tangent space of Xi in the 
random points. In the appendix we also show that these tangent spaces can be calculated 



EXPERIMENTAL RESULTS FOR THE POINCARE CENTER PROBLEM 



3 



using Frommer's algorithm. By using the inequality 

codim(Tx.:r) < codimX' 

where X' is a component of X passing through x we can eliminate all points on components 
of codimension at most c when estimating the number of components in codimension c + 1. 
This allows us to estimate also the number of components of larger codimension. 

The method of investigating many random points has minimal memory requirements and 
can be parallelized with almost no overhead. It also has the advantage of giving continuous 
intermediate results while running and can therefore be stopped and restarted with no loss. 

This characteristic p method differs from the one applied by Fronville |Fro97j . She reduces 
explicit focal polynomials modulo p to make factorisation and Grobner basis computations 
faster. For general cubic centers the explicit polynomials of high degree are unknown so 
Grobner basis and factorisation methods can not be applied. Even if the polynomials where 
known, they are to large to be factored with current methods, even over a finite field. 

We start by applying our method to homogeneous differential equations of degree 2 and 3 
where we recover the results known from literature. For general degree 3 differential equations 
we obtain 

Conjecture. Let X^ C A 14 be the algebraic set of inhomogeneous degree 3 Poincare differ- 
ential equations over C where all focal values vanish. Then Xao has 

(1) 1 component of codimension 5, 

(2) 2 components of codimension 6, 

(3) 4 components of codimension 7, 

(4) at least one component of codimension 8, 

(5) unknown numbers of components of codimension > 9. 

In Section 0] we investigate the known families of Darboux type Poincare differential equa- 
tions. Using a theorem of Christopher we exhibit Darboux type families whose closures form 
components of X^ with codimensions 5, 6, 7, 7 and 7. These are those with algebraic integral 
curves of degrees (4), (3, 1), (2, 2), (2, 1, 1) and (1, 1, 1, 1) respectively. 

In Section|S]we show that there exist two families of rational reversible systems whose closures 
form components of X^ of codimension 6 and 7. Namely these are the symmetric systems 
and Zoladek's family CR\\. Together with these results our experiments suggest that these 
components form a complete list up to codimension 7 and that Zoladek's conjecture holds up 
to this codimension. 

As a further check of our experiment we look a the other families of Zoladek's list in Section 
Eland show that all of them either lie on one of the above mentioned components or on com- 
ponents of codimension at least 8. For this we find new invariant conies for Zoladek's families 
Ci?5, Ci?7, CR\2 and CR±q and show that general members of these families are also of 
Darboux type. 

To conclude we show in Section [7] how our methods can be used for a heuristic estimate 
of the number of equations needed to define X^. This leads to the correct estimates for 
homogeneous equations of degree 2 and 3. For general degree 3 equations we do not yet have 
enough data to give a reliable estimate. Our experiments so far give a lower bound m(3) > 10, 
which doesn't quite reach the known bound m(3) > 11. To obtain heuristic evidence for a 
possible upper bound m(3) < 11 we would need to speed up our calculations by a factor of 
10000. 
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1. Preliminaries 



In this paper we write the differential equation 

, = P{x,y) 
Q{x,y) 

as P(x,y)dx + Q(x,y)dy = 0. If P and Q are polynomials of degree at most d we can 
homogenize P and Q with respect to a third variable z. Conversely we obtain a differential 
equation 

, = P(x,y,l) 
Q(x,y,l) 

from every homogeneous polynomial differential form P(x,y, z)dx + Q(x,y, z)dy. If not stated 
differently we will use the homogeneous formulation. 

Notation 1.1. Furthermore we denote by 

K an commutative ring (usually a field or Z), 

P 2 := ¥ 2 K the projective plane over this field, 

K[x,y,z] the coordinate ring of P 2 , 

A 2 C P 2 the affine plane where z / 0, 

P^o C P 2 the line at infinity where z = 0, 

K[x, y, z]d the vector space of homogeneous degree d polynomials, 

P,Q £ K[x, y, z]d two such polynomials, 

Pdx + Qdy the corresponding differential form on P 2 , 

V d = K[x, y, z]d © K[x, y, z]d the vector space of all such differential forms, 

{ptfyi z d ~ % ~i dx^x % yi z d ~ % ~i dy\i+j<_d the monomial basis of this vector space, 

PijiQij £ K the coordinates of V d with respect to this basis. 

Also we denote for any polynomial F £ K[x,y,z] the (formal) partial derivatives by F x , F y 
and F z . 

The differential equations considered by Poincare have a special form: 

Definition 1.2. Let Pdx + Qdy £ Vd be a polynomial differential form of degree d on P 2 . 
A point a £ P 2 is called a critical point of Pdx + Qdy, if P{a) = Q(a) = 0. Pdx + Qdy is 
called a Poincare differential form, if a = (0 : : 1) is a critical point, P x {a) = Q y (a) = 1 and 
Py{°) = Qx{o) = 0. We denote the vector space of all Poincare differential forms by Vp oincarfi . 
^Poincare ^ s a codimension 6 linear subspace of V d . 

Remark 1.3. The Poincare differential forms Pdx + Qdy G ^p incare correspond to the 
differential equations 

, = x + p(x,y) 
y + q(x,y) 

of the Poincare center problem, since the conditions above imply that P and Q are of the 
form 

P = xz d ~ x + terms with fewer z's and Q = yz d ~ l + terms with fewer z's 
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Definition 1.4. Let Pdx + Qdy be a Poincare differential form of degree d over a field of 
characteristic 0. One can then use Frommer's algorithm IA. II in the appendix to find a formal 
power series F E K[[x, y]] with 

\P(x,y,l) Q{x,y,l)J ^ Jy 

In this situation Sj(P,Q) is called the jth focal value of Pdx + Qdy. By Corollary IA.3I in 
the appendix Frommer's algorithm also implies that sj is polynomial on Vp oincar& and has 
rational coefficients. We call Sj E Q[pij,qij] the jth focal polynomial. 

Example 1.5. For Poincare differential equations of degree 3 the first focal polynomial is: 

Sl = §P02902 + 3P02PII - 5902911 + 3PIIP2O - 5911920 - §P20920 ~ P03 + ^912 ~ 3P2I + 930 

The size of the focal polynomials grows very fast. With current computer algebra systems 
one can not calculate Sk with k > 6 for degree 3 Poincare differential equations. For this 
reason we do not use explicit polynomials in this article. 

Definition 1.6. For E Q[pij,qij] we denote by 5^ the smallest common denominator of 
all coefficients, i.e. d^s^ E Z[py,gy]. We call Ij = (5xs\, . . . 5jSj) the i-th center ideal and 
-^oo = ($iSi)ieN the totai center ideal. 

The vanishing sets of these ideals are the i-th center variety X{ = V(I{) E ^p oincari5 and the 
total center variety = V(Ioo) C ^p oincari5 - Since Ii,Ioo C Z\pij,qij] these varieties are 
defined over every commutative ring K. 

Remark 1.7. By Corollary IA.5l in the appendix on can use Frommer's algorithm over a finite 
field to check 5jSj(P,Q) = without knowing SjSj explicitly. By Remark I A. 71 one can use 
Frommer's algorithm over K[e]/(e 2 ) to determine the tangent space to Xj in a given point 
Pdx + Qdy E Xj. For this it is again not necessary to know 5jSj explicitly. 



2. Counting Points 



In this section we explain how one can obtain heuristic information about a variety X C A n 
by evaluation its defining equations at random points. For an extended discussion about this 
method see |vBS04| . 

Definition 2.1. Let X C A n (¥ p ) be an algebraic variety. Denote the number of F p -rational 
points of X by \X(¥ p )\. Then 

p |A"(F P )| 

is called the fraction of¥ p -rational points of X in A n . 

Remark 2.2. If X has r irreducible components of codimension c and all other irreducible 

components have larger codimension then the Weil-Conjectures imply that 

/l\c 1 
7„(X) = r( — ) + higher order terms in — 
\pJ p 

We will estimate J P (X) statistically by evaluating the equations defining AT in a number of 
randomly chosen points. 

Definition 2.3. Let X C A n (F p ) be an algebraic variety. For a sequence S = (x%, . . . ,xn) 
of Fp-rational points in A ra (F p ) we call 

\{i\xi€X}\ 



%(X,S) 



N 
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the empirical fraction of¥ p -rational points. 

Remark 2.4. The distribution of ; y p (X,S) on the set of all sequences S of length N is 
binomial with mean p(j p (X, S)) = 7 P (A) and standard deviation 



This allows us to obtain an estimate of 7 P (A) and then of r and c by evaluating the equations 
of X in many random points. More information is obtained, if we also calculate the tangent 
space of X in these random points: 

Remark 2.5. Let X' C X C A n be an irreducible component, x € X' a point and Tx' x the 
tangent space of X' in x. Then 

codimX' > codimTx',x 

with equality for general points if X' is reduced. We therefore consider only points with 
codimTx',a; = c in estimating the number of components of codimension c. By the inequality 
above we disregard all points on components of codimension greater then c. 

These arguments lead us to 

Heuristic 2.6. Evaluate the equations of X in N random points Xj over F p and calculate 
the tangent spaces Tx, Xi m these points. Then estimate 

u r a- , x ^ I codimry.xj = c} c 
#{codim c components) sa — p 



with an estimated error 



J#{i | codimTx,^ = c} 

a iv ^ P 



c 



In this paper we have used a = 2. 

Caution 2.7. Let X c be the subvariety of X whose points have a tangent space of codimen- 
sion c. Then above heuristic means that statistically the hypothesis r ) p {X c ) = r(l/p) c can 
not be rejected with confidence of more than 4.6%. Algebraically this proves nothing, but 
gives a way to arrive at a reasonable conjecture about X. 

Remark 2.8. One might propose to estimate the higher coefficients in the power series of the 
Weil formulas using the same methods as above. Notice that in this case the error scales with 
(l/p) Cmin where c m i n is the codimension of the largest component. Using the infinitesimal 
information we have the better scaling (l/p) c . In our Experiment 13.31 the second method is 
23 2 ~ 500 times faster when estimating the number of codimension 7 components. 



3. Experiments 



To show that our heuristic gives useful results we applied it to the Poincare center problem 
with d = 2, d = 3 with p^j = q^j = for i + j = 2 (homogeneous case) and d = 3 without 
restrictions. For the first two cases the decomposition of X^ into irreducible components is 
known and agrees with the estimates of our heuristic. In the general d = 3 case our computer 
power is enough to estimate the number of components with codimension at most 7. 

Experiment 3.1. Over Fn we calculated the first 4 focal values of 100, 000 differential equa- 
tions with homogeneous nonlinearities of degree 2 and random coefficients using Frommer's 
Algorithm as described in the appendix. If all focal values vanished, we also calculated the 
codimension of the tangent space of X-j in this point. The results are collected in Figure ^ 
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Homogeneous Poincare differential equations 



Degree 2 

(p=ll, 100.000 points) 




2 3 
codimension of tangent space 



Degree 3 

(p=17, 68 million points) 




3 4 5 

codimension of tangent space 



Figure 1 . Our heuristic gives the correct results for homogeneous differential 
equations. 



Since we found not enough points with a tangent space of codimension or 1 to justify a 
component of this codimension we show them together with the codimension 2 points. 

Notice that our results agree with the well known fact that has 4 components in this 
case: The Hamiltonian differential equations (codim 2), the symmetric differential equations 
(codim 2), the component of differential equations with three integral lines (codim 2) and 
the component with an integral conic and and integral cubic in special position (codim 3) 
[Sch93] . 

Experiment 3.2. Over F17 we calculated the first 7 focal values of 68 million differential 
equations with homogeneous nonlinearities of degree 3 and random coefficients. If all focal 
values vanished, we also calculated the codimension of the tangent space of Xj in this point. 
The results are collected in Figure ^ Again we have shown the points with tangent spaces 
of codimension less then 3 together with those of codimension 3. 

Malkin Mal64 has first given center conditions in this case. More explicitly Lunkevich and 
Sibirskii have shown in [EB65 that 1^ decomposes into two ideals which they call condition 
I and II. Condition I defines a codimension 3 variety which on closer inspection decomposes 
further into the variety of Hamiltonian differential equations and the variety of symmetric 
differential equations. Condition II defines an irreducible variety of codimension 5. This 
agrees with our experiment. 

Experiment 3.3. Over ¥ p we calculated the first (p — 3)/2 focal values of several billion 
differential equations with general nonlinearities of degree at most 3 and random coefficients. 
If all focal values vanished, we also calculated the codimension of the tangent space of X^ p _ 3 y 2 
in this point. The results are collected in Figure [2j 

This experiment is consistent with the following conjecture 

Conjecture 3.4. Let X^ C A 14 be the algebraic set of inhomogeneous degree 3 Poincare 
differential equations over C where all focal values vanish. Then X^ has 

(1) 1 component of codimension 5, 
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General Poincare differential equations of degree 3 



p=17, 7 equations, 80 billion points 
p=19, 8 equations, 90 billion points 
p=23, 10 equations, 330 billion points 



I 



I I I I 

5 6 7 8 

Codimension of tangent space 

Figure 2. The higher result for p = 17 in codimension 7 is expected, since 
in this case we only consider Xj. The 8th focal value could be nonzero on a 
general point of one component. Similarly we expect a higher result for p = 19 
in codimension 8. In codimension 8 we have only 2 measurements, since one 
can not evaluate 5$s$ over F17 with Frommers algorithm. 



(2) 2 components of codimension 6, 

(3) 4 components of codimension 7, 

(4) at least one component of codimension 8, 

(5) unknown numbers of components of codimension > 9. 



In the following sections we will identify known families of differential forms with a center 
whose closure is a component of X^ with codimension at most 7. The number of such families 
agrees with our experimental results. This provides new evidence for Zoladek's conjecture up 
to codimensin 7. 



4. Darboux Centers 

Let Pdx + Qdy be a differential form of degree d with algebraic integral curves whose degrees 
sum up to d+1. If these integral curves are in general position it is well known that Pdx+Qdy 
admits an integrating factor KC93]. If Pdx + Qdy is a Poincare differential form this means 
that Pdx + Qdy has closed integral curves near the origin. 
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In this section we consider for each partition \i of d + 1 the variety of Poincare differential 
forms with algebraic integral curves of degrees [i in general position. For d = 3 we calculate 
their dimension and show that their closures form components of X^. 

Notation 4.1. For the purpose of this section we consider 

F € K [x, y, z] a homogeneous Polynomial of degree e 
Cp := {F = 0} the corresponding algebraic curve in P 2 

Remark 4.2. is an integral curve of Pdx + Qdy if and only if PF y — QF X = on Cp. 
This is the case if and only if there exists a unique homogeneous polynomial Kp 6 K[x, y, z] 
of degree d — 1 such that 



FFy QF X 



-FKp. 



We put the minus on the right side to simplify the formulas later on. 
Definition 4.3. In this situation Kp is called the cofactor of F. 

The importance of algebraic integral curves and their cofactors is given by 

Theorem 4.4 (Darboux). Let F±, . . . , F n € K[x, y, z] be homogeneous Polynomials, such that 
the corresponding algebraic curves C\, . . . , C n are integral curves of Pdx + Qdy with cofactors 
K\ , . . . , K n . If there exist scalars an £ K such that 



Qx Py — ^ ^ C^iF-i 



i=l 



then 



is an integrating factor for Pdx + Qdy . Such differential forms are called Darboux integr able. 

Proof. |Dar78| □ 

If the algebraic integral curves have the right degrees, one doesn't even have to consider the 
cofactors: 

Theorem 4.5 (Christopher, Zoladek). Let C\, . . . ,C n C P 2 be smooth algebraic curves of 
degrees e±, . . . , e n with 

(1) no curve tangential to the line at infinity, 

(2) no two curves tangential to each other, 

(3) no two curves meeting the line at infinity at the same point, 

(4) no three curves meeting in the same point. 

Then the vector space of degree d = e\ H \- e n — 1 differential forms Pdx + Qdy that have 

these integral curves is n dimensional. Further more all these differential forms are Darboux 
integrable. 

Proof. Pdx + Qdy has C\ , . . . , C n as integral curves with cofactors K\ , . . . , K n if and only if 
(Ki, . . . , K n , -Q, P) • M = where 





fF 1 . 


. \ 


M = 


. 


• F n 




Fix ■ 


TP 

± nx 




\F ly . 


■ Fny/ 



10 



V. BOTHMER 



and the Fi are defining equations of the Cj. The conditions of the theorem are satisfied if 
and only if M has full rank everywhere on P 2 . In this case a presentation of the kernel of 
M is given by the Buchsbaum- Rim-complex Eis95( Apendix A2.6]. This complex gives the 
number of degree d differential forms and writes them in terms of n x n minors of M . A short 
calculation then gives that Q x — P y is in the span of K\ , • • • , K n . For a proof involving only 
Hilbert's Nullstellensatz and no complexes see |CCZ02| Theorem 1] and KC93, Theorem 
1]. " □ 

We now pass from general differential forms Pdx + Qdy to Poincare differential forms: 

Lemma 4.6. Let Pdx + Qdy a Poincare differential form that is Darboux integrable with 
respect to the algebraic integral curves C%, . . . , C n . Let F\, . . . , F n be the defining equations and 
K\, . . . , K n the cof actors. I/o£F 2 is a point not lying on any Cj satisfying P(a) = Q(a) = 
then Q x (a) = P y {a)- 

Proof For each curve Ci we have 

P{a)F iy {a) - Q{a)F ix + F^K^a) = => K^a) = 
since P(a) = Q(a) = and Fi(a) ^ 0. Consequently 

n 

Q x {a) - P y {a) = a t Ki{a) = 
i=0 

since Pdx + Qdy is Darboux integrable with respect to the Ci. □ 

Definition 4.7. We denote by Y\ with A = (Ai, . . . , A n ) h d+1 the variety of differential forms 
that have n integral curves of degree A, satisfying the conditions of Christopher's Theorem. 
For the subvariety of Poincare differential forms we use the notation X\ = Y\ fl Vp oincar £. 

Proposition 4.8. Let d = 3. Then 



(1) Xu\ is a codimension 5 component of X^, 

(2) A( 31 ) is a codimension 6 component of X^, 

(3) X (2,2)1 -^(2,1,1); ^(1,1,1,1) are codimension 7 components of Xoo. 

Proof. The variety of degree 4 smooth plane curves, that do not pass through the origin 
is 14 dimensional. By Christopher's Theorem this variety contains an Zariski open subset, 
such that for each curve of this subset there is a 1 dimensional family of differential forms 
Pdx + Qdy that have this integral curve. For degree reasons the cofactor in this case has 
to be 0. By Christopher's Theorem these differential forms are Darboux-integrable with 
Qx — Py = 0. This means that Pdx + Qdy is exact and there exists a polynomial H of degree 
4 with P = H x and Q = H y . Therefore each such differential form has a 1-dimensional family 
of algebraic degree 4 integral curves. 

Furthermore the Poincare differential forms Vp oincari s are defined by 6 linear equations in the 
space of all differential forms of degree 3. By Lemma l4.fil at most 5 of these are independent 
on the set of Darboux-integrable differential forms. Combining these dimension counts we 
obtain 

dimX (4) > 14 + 1-1-5 = 9. 
Since dimVp oincar ^ = 14 this means codimX( 4 ) < 5. On the other hand 

(x-9x 2 +2xy+2y 2 +2xy 2 -y 3 )dx+(y+x 2 +4xy-6y 2 -2x 2 y-3xy 2 )dy 

is a Hamiltonian Poincare differential form with potential 

(1 /2)x 2 +(l/2)y 2 -3x 3 +x 2 y+2xy 2 -2y 3 +x 2 y 2 -xy 3 . 
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It has a 1-dimensional family of smooth plane quartic integral curves that do not pass through 
the origin and the tangent space of X\q at this point over F23 is 5. By semicontinuity 
arguments this proves that codim Xu) > 5 and that Xu\ is a component of X^, . 

The variety of smooth degree 1 and degree 3 curves have dimension 2 and 9 respectively. The 
variety of pairs that intersect transversally is then 11 dimensional. For an open subset of 
these pairs there is a 2 dimensional set of differential forms that have these integral curves. 
Again by Lemma 14.61 the codimension of the variety of Poincare differential forms in this set 
is at most 5. Since 

(x—10x 2 +2xy—2y 2 —8x 2 y—llxy 2 +3y 3 )dx+(y—5x 2 —4xy—3y 2 —x 3 —9x 2 y+10xy 2 —8y 3 )dy 

is a Poincare differential form with exactly one degree 1 and one degree 3 integral curve 

10y+z=0 and 10x 3 -llx 2 y+6xy 2 -9y 3 +10x 2 z+2y 2 z+2yz 2 +z 3 =0 

over i*23 satisfying the conditions of Christopher's Theorem, we have 

dimX (3jl) > 11 + 2- 0- 5 = 8. 

Since the tangent space of X\q at this point has codimension 6, this proves that ^(3,1) is a 
codim 6 component of X^. 

For the remaining components we observe that the variety of degree 2 plane curves is 5 
dimensional and that 

7 = (5 + 5) + 2- 5 
= (5 + 2 + 2) + 3- 5 
= (2 + 2 + 2 + 2) + 4-5. 

Now the Poincare differential form 

{x-2x 2 -llxy+9y 2 +3x 3 +2x 2 y-llxy 2 +lly 3 )dx+(y+4x 2 -3xy+y 2 +2x 3 -llx 2 y+llxy 2 -7y 3 )dy 

has two integral conies 

-6x 2 -7xy-10y 2 +6xz-llyz+z 2 =0 and 2x 2 +8xy+lly 2 ~4xz-8yz+z 2 =0 

in satisfying the conditions of Christopher's Theorem over F23. The Poincare differential 
form 

(x-3x 2 +llxy+5y 2 +8x 3 +2x 2 y+8xy 2 —4y 3 )dx+(y-5x 2 -3xy-lly 2 +7x 3 +5x 2 y-3xy 2 )dy 

has two integral lines and one integral conic 

-6x+z=0, -3x+3y+z=0 and x 2 +5xy-2y 2 -4xz-5yz+z 2 =0 

in satisfying the conditions of Christopher's Theorem over F23. The Poincare differential 
form 

(x-llx 2 +3xy-2y 2 -3x 3 +9x 2 y-5xy 2 -lly 3 )dx+(y+5x 2 +8xy-lly 2 +4x 3 -4x 2 y+8xy 2 +7y 3 )dy 

has four integral lines 

3x+8y+z=0, -2x+9y+z=0, x+5y+z=0 and 10x-6y+z=0 

in satisfying the conditions of Christopher's Theorem over F23. 

Moreover the tangent spaces of X±q in these points over F23 are all of codimension 7. □ 
Remark 4.9. All examples in this section were found with the program |vBC 05 . 
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5. Rationally Reversible Systems 
A second type of centers has been considered by Zoladek : 

Definition 5.1. A differential form Pdx + Qdy is called rationally reversible, if there exist 
a rational map 

$ . c 2 -> C 2 

and a second differential form P'dx + Q'dy such that 

<S>*(P'dx + Q'dy) = n(Pdx + Qdy) 
with // a suitable polynomial. 

This definition is useful because of 

Proposition 5.2 (Zoladek). Let fi(Pdx + Qdy) = (P' dx + Q' dy) be a rationally reversible 
differential form over R and O G M? a ramification point of . Furthermore doenote by T$ 
the ramification curve of & and <&(T$) the branch curve. If 

(1) there exists a neighborhood U C M 2 of O such that the boundary of $(U) contains 
part of the branch curve 3>(r$), 

(2) P'dx + Q'dy does not vanish at &(0) and 

(3) the integral curve of P'dx + Q'dy through &(0) is tangent to the branch curve 3>(r$) 
from the outside o/ <&(£/), 

then O is a center of Pdx + Qdy. 

Proof. (Zoladek) In this case the integral curves of Pdx + Qdy close to O are preimages of 
compact curves in &(U). □ 

Definition 5.3. In the situation above, let Tf id C T$ be the union of those components of 
r$ that contain O. Tf id is then called the fold curve of <I>. 

The first examples of rationally reversible differential forms are those which are symmetric 
with respect to the x axis. These were already considered by Poincare. In this case $ = 
(x,y 2 ). In Zoladek's list this ist the family CR\. 

We now show, that the Poincare differential forms which are symmetric to a line though the 
origin, form a codimension 6 component X sym of X^. 

Proposition 5.4. Let X x _ ax i St c C ^p i ncare - be the variety of complex valued Poincare dif- 
ferential forms that are mirror symmetric with respect to the x-axis, i.e. those Pdx + Qdy 
that satisfy P(x, —y, z) = P(x, y, z) and Q(x, —y, z) = —Q(x, y, z). Then X x ^ ax i s C X 



Proof. X x _ ax i s ^c is a vector subspace of ^p ; ncar ^, since symmetry with respect to the x axis 
is a linear condition. For real symmetric differential forms X x _ ax i s ^ C X x ^ ax i S) c the focal 
values Sj vanish, (see for example Pre03, p. 42, Satz 5.2.1]). Since the Sj are holomorphic 
on X x — ax i s ^ they have to vanish on all of X x — ax i s 

,c- □ 

Remark 5.5. The same result can also be obtained analytically by introducing complex 
polar coordinates |vW05| . 

Remark 5.6. Over the real field mirror symmetry implies stability even for differential 
equations of the form 

x 2n ^+p(x,y) 

"' y 2n ~ l + q(x, y) 

|Pref)31 p. 42, Satz 5.2.1]. 
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Proposition 5.7. Let X sym C C Vp oincarf s be the variety of mirror symmetric Poincare 
differential forms of degree 3. Then X sym is a component of X^ and cod\mX sym = 6. 

Proof. X x - ax i s is defined by the linear conditions pij = for j odd and qij = for j even. 
Therefore codim X x _ ax i s = 7. 

By rotating the symmetry axis about the origin, i.e. by applying an element of SO(2,C), we 
obtain a 1-parameter family of mirror symmetric differential forms in for each element 
Pdx + Qdy e X x 

—axis- Since for a general Pdx + Qdy (E X x — ax i s this family is not constant, 
this shows codim X sym < 6. 

Now 

(x + llx 2 + y 2 —9xy 2 )dx + (y—Qxy + x 2 y — 5y 3 )dy 

is mirror symmetric. In characteristic 23 it can be checked with Prommer's algorithm that 
codim T X) x 10 = 6. Since X sym C X\q, this shows codim X sym > 6. It also shows that the 
component of Xoo that passes through x has codimension at least 6. Since X sym is irreducible 
of codimension at most 6 and contains x, it must be this component. 

□ 

Next we consider Zoladek's family CR\\. 

Definition 5.8. Let $ = (A 2 /B, A/C) with A,B and C linear polynomials in x and y, and 

P'dx + Q'dy = —2x(ky + Ix + mxy)dx + y(ky + nx + py 2 + qxy + rxy 2 )dy. 

Then 

Ycmi = {Pdx + Qdy\3A, B,C, k,l,m,n,p,q,r, fi: fi(Pdx + Qdy) = <$>* (P' dx + Q' dy)} 
is called Zoladek's family CR\\. We also set Xcmi ■= Ycmi D ^p mcare- 



Proposition 5.9. Xcmi is a codim 7 component of X^. 

Proof. Since Ycru is rational and the focal values are holomorphic functions on the open 
set of differential forms with an elementray center that vanish on the real points of Ycru by 
Zoladek's Theorem 1 in |ZoI94 they must also vanish on all of Ycru- 



Also Yqru is a codimension 7 family of differential forms by the same theorem. Therefore 
^Ci?n has codimension at most 7. Now 

(— 8x 3 +12x 2 y+6xy 2 — 18y 3 +x 2 — 2xy— 6y 2 +x)dx+(5x 3 +3x 2 y— Uxy 2 +lly 3 +2x 2 — 15xy— 15y 2 +y)dy 

is an element of Xcru with a codim 7 tangent space in X\j over F37. This proves the 
proposition. □ 

Remark 5.10. For a random subset of the codim 7 points we found in our experiments we 
also checked wether they lie on Xcru- This was the case for about one quarter of the points, 
which is consistent with our conjecture that there are exactly 4 codimension 7 components 
of which Xcru is one. 

6. Zoladek's List 

In an amazing work Zolade k classifi ed all rationally reversible cubic systems t hat are not 



Darboux integrable ( Zol94 , Zol9£|). A list of the 17 such centers is presented in |ZoI94j and 
some minor mistakes are corrected in |Zol96 . In this section we compare these results to our 
experiments. 
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Proposition 6.1. Let Ycm C V s be the family CRi from Zoladek's list, and Xcm = Ycm H 

Vpoincare ^ e SU0S P ace of Poincare differential forms in this family. Then 

• codimXcTjj > 6 for i = 1, 

• codim Xcm > 7 for i = 5,7, 11, 12, 16, 

• codim > 8 /or % = 2, 4, 6, 8, 13, 14, 

• codim Xcm > 9 /or i = 3, 9, 10, 15, 

• codim Xcmi > 10 /or i = 17. 

Proof. Points on Xcm that have the have tangent spaces of the above codimension can be 
found by the following method. Choose a random differential Form in a; € Ycm ° ver a finite 
field F p . If there is a F p -rational point on the fold curve, where the corresponding differential 
form degenerates, apply a change of coordinates 4> that moves this point to the origin. If 
the linear part of x' = 4>(x) admits a F p -rational coordinate change ip such that ip(x') is 
in Vp oincar( s set y = ip(x') and calculate the tangent space with Frommer's algorithm. This 
methods usually gives an example after 4 or 5 trials. The proposition then follows from 
semi-continuity. □ 



Remark 6.2. In Zoladek's list of Zol96| there are some minor misprints. For CR% one has 
to replace (2a 2 — b) in the formula for F by (2a — b 2 ). For CR5 the fromulas for x and y 
have to be exchanged and —I has to be replaced by —ly . For CR12 the ter m —pqT must be 



replaced by —yqT. In every case either the list in [Zol94| . the one in [Zol96 or a combination 
of the two is correct. 



Remark 6.3. Because of the new formula for r\ given in Zol96j the fold curve in this case 
changes to T = {x — 2ay + 1 — a = 0}. 

Remark 6.4. Because codim(Xcm C Vp i n care) — co ^ m 0^cm C V 3 ) our examples for the 
family CR\§ show that the codimension 5 given in Zoladek's list can not be correct. 

We will now take a closer look at those families that could be of codimension 7 in the space 
of Poincare differential forms. As in Zoladek's papers we set 

T = x + y + c and T2 = ax 2 + bxy + cy 2 + dx + ey + 1. 

We have found some previously unknown conic integeral curves in some of Zoladek's rationally 
reversible families, proving that they are also Darboux integrable: 

Proposition 6.5. A general differential form of Zoladek's family CR\q 

Pdx+Qdy=— (^—y(ky+lx) — (py+qx)*T2 — (2ax+by+d)(—qx 2 + (n—p)xy+rny 2 )j dx 
+ (^—x{ky+lx) — {my+nx)T2+(bx+2cy+e)(—qx 2 + {n—p)xy J rmy 2 )} dy 

has a conic integral curve 

{(mg — np)T2 + (kq — lp)x — (kn — lm)y = 0} 



and is Darboux integrable with Xcriq C -X"(2,i,i) • 

Proof. That the conic above is an integra l curve can be checked by a straight forward calcu- 



lation. By Zoladek's Remark 1 in |Zol94j the differential forms of this system also have two 
integral lines given by the equation qx 2 + (p — n)xy — my 2 = 0. By substituting random values 
of F p for the parameters of this family one easily obtains examples that satisfy the conditions 
Christophers Theorem 14.51 over F p . The proposition then follows by semi-continuity. □ 



EXPERIMENTAL RESULTS FOR THE POINCARE CENTER PROBLEM 



15 



Proposition 6.6. A general differential form of Zoladek's family CR\ 2 

Pdx+Qdy=-( y -ky 2 +2{q-m)xy+2rxT-2ly' i -yqT-rT 2 -2ny 2 T-2pyT 2 ")dx 
+2{jjn-k-q)xy-{m+r)xT+ly' A +(p-n)yT 2 +(n~l)y 2 T-pT 3 )dy 

has a conic integral curve 

{a 2 ox 2 + anxy + a 02 y 2 + a 1Q x + a 01 y + a 0Q = 0} 
where the a%j satisfy linear system of the equations 

( -c 

— 2pc — m — r 1 1 It 

—2nc — k + 2m — q + r —m + l/2q — r 

—2lc+k — m + q —l/2k + m — q —m + q I 






r 





P 

n 



1 0\ 









V 















V 



Furthermore Pdx + Qdy is Darboux integrable and Xcru C Xr 2 . 



( O2o\ 

an 
ao2 
aio 
aoi 
\ a oo/ 



0. 



Proof. Same as for Proposition 16.51 with Zoladek's integral lines given by the equation ky 2 + 
qTy + rT 2 = 0. (In Zoladek's paper the equation is ky 2 +pTy-\-qT 2 = 0, but this is a misprint 
as can be checked by a straight forward computation). □ 

For the families CR-j and CR§ the integral curves we found do not satisfy the conditions of 
Christophers Theorem. To prove that they are subfamilies of -^(2,2) we need the following 
technical lemma: 

Lemma 6.7. Let Pdx + Qdy be a differential form of degree 2>, 

(Fx Gx 
Fy Gy 
aF bG, 

a matrix with F, G quadratic polynomials and a, b E C\{0} such that 

(1) (-Q, P, Q x -P y )-M = and 

(2) M drops rank in only finitely many points o/P^. 



Then Pdx + Qdy is Darboux integrable and an element of Xm2) ■ 

Proof. In this situation F and G are integral curves of Pdx + Qdy with cofactors a(Q x — P y ) 
and b{Q x — Py). In particular Pdx + Qdy is Darboux integrable. 

If M drops rank in codimension 2 then the syzygies of M are generated by 







Gy 




Fx 


G x 




F x 


Gy 


( 


aF 


bG 


1 


aF 


bG 




Fy 


Gy 



For degree reasons the syzygy (—Q,P,Q a 
Q = -c 



F O 

JTy yjy 

aF bG 



P y ) has to be a scalar multiple of s, in particular 

F x G x 
C aF bG ' 



and 



P 



We now deform F and G to F and G such that C p and Cq satisfy the conditions of Christo- 
pher's Theorem. This is possible, since these conditions are Zariski open on the set of all 
pairs of conies. Then Cp and Cq are integral curves of Pdx + Qdy with 



Q 



F G 

r y ^y 

aF bG 



and 



P 



F x G x 
aF bG 
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and cofactors 



Ft, G, 



F x G, 



F G 



ca ~ ~" and cb 

ty Gy 

By Christophers Theorem Pdx+Qdy G X(2,2) is Darboux integrable. Now Pdx+Qdy € -^(2,2) 



deforms to Pdx + Qdy, which is therefore in -X"(2,2)- 

Proposition 6.8. A general differential form of Zoladek's family CR7 

Pdx+Qdy=— (nx+ky+nT+(m—l)x 2 y+px 2 T+mxyT+pxT 2 ^ dx+x(—(n+k)+(l—m)xy—(l+p)xTj dy 

has a conic integral curve 

Cf = { (lp — mp)x 2 + (Im — m 2 + lp — mp)xy + c(lp — mp)x + (—mn + kp) = 0} . 



□ 



and is Darboux integrable with Xcrj C -^(2,2) • 

Proof. That the conic Cp above is an integra l curve can be checked by a straight forward 
calculation. By Zoladek's Remark 1 in Zol94 the differential forms of this system also have 



another integral conic Cg with G = k + ITx. With a 



l—m 



and b 



we obtain a 



21— m 21— m 

matrix M as in Lemma 16.71 Condition (1) of Lemma 16 . 71 follows from a direct computation 
and condition (2) can be verified for random values of the coefficients over a finite field, since 
the rank condition is Zariski open on the variety of all such matrices. The claim then follows 
from Lemma 16.71 □ 

The most involved argument for Darboux integrability is needed for the family CR5: 
Proposition 6.9. A general differential form of Zoladek's family CR5 

Pdx+Qdy=(^—kxy 2 —ly+mx 2 y—(nxy+p+qTx)(2x+y+c)^ dx—x(l+p+mcx+{k+n)xy+mx 2 +qTx)dy 

has two conic integral curves 

{a2ox 2 + a\\xy + ca 2 ox + a o = 0} 
where a^- are such that the matrix 



M 



( ^20 Q ~m -q\ 

—an m — n — q k + n + i 

a 00 P -l-p 

\ -a 2 o an / 

has rank at most 2. Pdx + Qdy is Darboux integrable with respect to these quadrics and 

XcFCo C -^"(2,2) ■ 

Proof. Indeed for a general member of Zoladek's family CR5 we may after row and column 
operations assume that M has the form 

/0 1 0\ 
1 
LOO 
\Q 0/ 

with L linear in a^- and Q quadratic in a^. So M drops rank in two points and we have 
two quadrics. With a computer algebra program one can check that the equations given 
define an ideal, that contains the ideal that describes the set of all quadric integral curves. 
Consequently the quadric of the proposition are integral. 
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Now consider the matrix 

(2ym — yn 
2xm — xn 
—xqtc + xym — xyn — x 2 



2mtc — ntc + Axm — 2xn 


xyq — pt 2 xmtc + xqtc + xyk + x 2 m + xyn + x 2 q + xyq + It 2 + pt 2 



which we obtained as syzygy matrix of (— Q, P, Q x — P y ). It does indeed satisfy (— Q, P, Q x — 
Py) ■ M = 0. To show that it can be transformed into the form of Lemma lH. 71 we multiply the 
last line by a and subtract its derivatives from the first two lines. Looking at the coefficients 
of the resulting linear entries we obtain a 2-matrix 



M 





2q 




Anna - 


- 2na — 


2m - 


2q 


2ma — 


net — m - 


- n + q 




—k — n 


- 9 




2ma — 


net — m - 


-n + q 




~k — n 


- <? 






















qc 




2mac 


— nac - 


- mc — 


qc 



















This matrix has a kernel if and only if the first 2x2 minor 

<j) = (8m 2 — 8mn + 2n 2 )(a — a 2 ) — 2m 2 + 2mn — 2kq 

vanishes, since all other rows are dependent on the first two. For generic choices of k, m, n 
and q the polynomial 4> has two zeros a\ ^ OL2 ^ 0. For each of these values M' has a kernel 
Pi = (Pn,f3i2) T - By construction we obtain 



M 



Pn 
P12 



P21 
P22 



F„ 



-F 




with appropriate quadric polynomials F and G. 



By substituting random values for c, k, m, n and q over a finite field ¥ p we easily find 
an example where a\ and «2 are F p -rational and the constructed matrix satisfies the rank 
condition of Lemma 16.71 The proposition then follows by semi continuity. □ 

Corollary 6.10. The only family of rationally reversible cubic centers of codimension 7 in 
^Poincare whose general member is not Darboux integrable with integral curves of degree at 
most two is CR\\. 



Proof. By looking at random points of CR\\ over ¥ p one can easily find an example that 
is not Darboux integrable with integral curves of degree at most two. CR\ is the codim 
6 component of symmetric differential forms. All other families of Zoladek's list are either 
Darboux integrable or have codimension at least 8. □ 

Remark 6.11. It would be interesting to find a geometric explanation for the extra conic 
integral curves we have found. 



7. Necessary equations 

As a last application of our methods we estimate the number of focal polynomials needed to 
define X^. 

Proposition 7.1. Let x € A n be a point with si(x) = ■ ■ ■ = Sfc(x) = mod p, Sfc+i(x) 7^ 
and codim Tx fci a; = k where = V(s±, . . . ,Sk)-Then there exists a component of X^ that is 
not contained in Afc +1 in characteristic 0. 
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Proof. Since the maximal codimension of is k, x lies on a component X' of X^ in char- 
acteristic p that is smooth in a neighborhood of x and of codimension k. By a theorem 
of Schreyer ( |Sch96j . vBEL05]), this component lifts to a component X' over an algebraic 
number field. Now x does not lie in X^+i considered as a scheme over Z. Therefore Sk can 
not vanish on all of X' C □ 

Example 7.2. For 

Pdx+Qdy=(x-4x 2 +7xy+lly 2 +7x 3 +7a; 2 y-6xy 2 -3y 3 )rfx+(y-3x 2 -5xy-% 2 -10a' 3 ~4a 2 y+llxy 2 -6y 3 )dj/ 

we have si(P,Q) = ■■■ = sg(P,Q) = mod 23 and siq(P,Q) ^ mod 23. Therefore 
771(3) > 10. This example was found with the program |vBC05j . Notice that c urrently the 
best known bound is m(3) > 11 which requieres a more ingeneous proof Zol95 , |Chr 05b . 

For upper bounds we have to fall back on heuristic methods: 

Heuristic 7.3. If X^ ^ X^+i there exists a component X' C X^ that is not completely 
contained in X^i- Since all components of X^ have codimension at most k one would 
expect to find at least one point of X' with probability 

l_(l_l XiV 

V p k 

If this probability is large and we do not find such a point after N trials we conjecture 
Xk = Xk + i. 

Example 7.4. In Experiment 13.11 we expect to find at least one point on A3 — X4 with 
probability 

x 100,000 



1 - (^1 - J ~ 99, 945%. 

if such points exist. Since we didn't find any points on X3 — X^ this leads us to estimate 
m(2) = 3. Indeed this is the correct value. 

Example 7.5. In Experiment 13.21 we would expect to find at least one point on X5 — Xq 
with probability 



j x 68,000,000 

1-I1-— =•) « (100 - 1.5 x 10" 19 ) 



70. 



17 5 

if such points exist. Since we did not find any points on X5 — Xq we conjecture that 
77i<(3homogeneous) = 5. Indeed the intersection of condition I and condition II in |LS65j is 
generated by the first 5 focal polynomials. 

Example 7.6. To obtain heuristic evidence for m(3) = 11 one would need to evaluate 

N = ln( °-° 5 1 ) - « -ln(0.05)29 n « 3.7 x 10 16 
ln(l - 2gir) 

random points. Unfortunately this about 10000 times more than we can currently manage. 



Appendix A. Frommer's algorithm 



Hans-Christian Graf v. Bothmer and Martin Cremer 

In this appendix we state Frommer's algorithm and show how one can use it to define the 
focal polynomials SjSj, evaluate them over finite fields and calculate tangent spaces to their 
vanishing sets without ever writing them down explicitly. 



EXPERIMENTAL RESULTS FOR THE POINCARE CENTER PROBLEM 



19 



Algorithm A.l (Frommer, Moritzen). Let Pdx + Qdy be a Poincare differential form on 
PjL For n < 2k + 2 calculate q n _/ for < / < n and a n _^ according to the following 
formulas 

Q, n _z := ^2 {-{ n ~ 1 ~ J + l)Pijai-i,n-i-j+i + + l)qijai- i+ i, n -i-j) 

2<i+j<n 

n—l—2 



Q>n-U '■— \ 



-E 

8=1 

a(min(^i,^i),n) 



n — I 

c 2i-l,n-(2i-l) If n — 2 j 



n — / 



n 



2j 







2 min(7, n — I) 

2 

E 



Z-l 
2 

E 



C n -2i,2i 



^i=0 

i-l i-3 
2 2 



n 



a(i, n) 
n - 2j - 1 



E 



g "-2j,2j 

a(j, ra) 



Z - LJ - 2j + l 

n — l — 1 n — l — 3 

C2i,n-2i TT n-2j -1 



E 



n 



n — I 

i=0 j=i 



2j + l 

n-(2k- 1) 
(2k - 1) 



for /, n even, / < n 
for I = n even 

for I odd and n even 

for I, n odd 

for I even and n odd 



0. Also set Ojj = cy 



if 



where the constants a(i, n) are defined as 

% 

a(i,n) := ]J 
fc=i 

As start values use ci2,o = 00,2 = l)Oo,o = °o,i = °i,0 = °i,i 
either i < or j < 0. Then 

fc-fT 

/ n 1 c 2fc+2-2i,2i 

i=0 v ' ; 

is called the k-th focal value of Pdx + Qdy. 

Proposition A. 2. If char K = Frommer's algorithm is well defined and the formal power 
series F(x,y) = Ylij a ij x% y^ ^ -^[[^y]] satisfies 



det 



Fx(x,y) F y (x,y) 
P(x,y,l) Q(x,y,l) 



E 

3=1 



Si(x 2 J+ 2 + y 2 J+2) 



Proof. In characteristic zero all denominators that appear in Frommer's algorithm are invert- 
ible. Furthermore the formula for q n _j involves only a« with i + j < n and the formula for 
a n -u involves only with i + j < n. Therefore Sk(P, Q) is well defined. Moritzen shows in 
MorOO that F satisfies the above formula. 

The idea and the first few steps of this algorithm appear in |Fro34j . The explicit formulas 
above we re fou nd by Moritzen IMorOOl . For our implementation of Frommer's algorithm we 
followed HOI . We used a random number generator from |PTVF92] . The source files of our 
program can be obtained from |vBC05| . □ 

Corollary A. 3. The function : V d — ► K is polynomial in pij and qij with rational coeffi- 
cients, i.e s k G Q\pij,qij]- 



Proof. In Frommer's algorithm all formulas are algebraic and only integral numbers occur as 
denominators. □ 
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Corollary A. 4. If char K = p then Sk is well defined if 2k + 2 < p. 



Proof. The denominators in the formulas for q n _; and a n ^ij and Sk in Frommer's algorithms 
are products of natural numbers less or equal to n < 2k + 2. □ 

Corollary A. 5. Let ¥ p be a finite field of characteristic p > 2k + 2 and Pdx + Qdy a 

differential form over Z. Zf s^p (P,Q) is the result of Frommer's algorithm over ¥ p and 
Sk,i$(P-Q) the result of Frommer's algorithm over Q, then 

s k ,F p (P, Q) = 5 k s ktQ (P,Q) = 

where 8k the smallest common denominator of the focal polynomial s k ,Q £ Q[Py')9ij] an d 
Pdx + Qdy is the reduction of Pdx + Qdy to ¥ p . 

Proof. By the argument above, 5k is not divisible by p. □ 

Remark A. 6. This allows us to find points on Xi over a finite field with Frommer's algorithm 
without knowing the explicit polynomials 5k s k- 

Remark A. 7. For a Poincare differential form Pdx+Qdy over a field K and vector (P', Q ; ) £ 

^Poincare We haVe 

Sfc (P + eP', Q + eQ') = s k (P, Q) + es' k (P, Q) G K[e)/(e 2 ) 

where s' k is the formal derivative of Sk in the direction (P',Q'). Since Frommer's algorithm 
works over any ring of characteristic p > 2k + 2 this allows us to calculate the tangent space 
of Xk in the point Pdx + Qdy without knowing the explicit polynomials SkSk- 

Remark A. 8. To speed up the search for rational points on X k over ¥ p , we used the first 
focal polynomial (see Example 11.5(1 to calculate q^o from the other values. This effectively 
multiplies the number of points checked by p. The numbers reported in the main part of 
this paper are effective numbers, not actual numbers. Alternatively we could have used the 
actual number of points and calculated the codimensions in the hypersurface s\ = and then 
added 1 to obtain the codimension in V^- n -. The results are the same. 
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